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Abstract 



We calculate the 2-loop partition function of QCD on the lattice, using the 
Wilson formulation for gluons and the overlap-Dirac operator for fermions. 
Direct by-products of our result are the 2-loop free energy and average pla- 
quette. 

Our calculation serves also as a prototype for further higher loop calcula- 
tions in the overlap formalism. 

We present our results as a function of a free parameter Mq entering the 
overlap action; the dependence on the number of colors N and fermionic 
flavors Nf is shown explicitly. 
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I. INTRODUCTION 



In recent years there has been great progress toward the implementation of chirahty 
preserving regularizations of gauge theories coupled to fermions on the lattice. In the two 
most widely studied approaches, overlap fermions and domain wall fermions (for recent 
reviews, see [1,2]), numerical simulations are already producing promising results. 

Alongside with numerical simulations, a number of perturbative calculations have also 
been carried out [3-6], which allow for a meaningful comparison between lattice results 
and continuum quantities. However, unlike the more conventional fermionic actions, the 
complexity of the overlap and domain wall actions has hampered the scope of perturbative 
results, and has kept the order of perturbative calculations down to 1 loop. 

In this paper we perform a first 2-loop calculation using the overlap-Dirac fermionic 
action. The quantity we compute is the expectation value of the plaqucttc, in SU{N) 
gauge theory coupled to Nf massless fermionic species. Aside from serving as a longtime 
testing ground for lattice perturbation theory, the plaquette operator has appeared in a 
variety of lattice investigations. One of its earliest uses has been in attempts to extract 
the quadratic gluonic condensate [7,8]. Another use regards the definition of an effective 
coupling constant [9,10], in an effort to improve the scaling properties of observables; this 
definition leads also to a modification of asymptotic scaling, which involves higher-loop 
coefficients of the plaquette. Finally, the plaquette is a standard ingredient in studies of the 
quark- ant iquark potential (see, e.g., [11] and references therein). 

In standard notation, we write the gluonic and fermionic parts of the action as follows: 



As usual, f3 = 2N/gQ is the lattice coupling and a is the lattice spacing; □ stands for the 
plaquette and Ux,x+n for the link variable joining neighboring sites in the jl direction: x and 
x + ji. The Neuberger-Dirac fermionic operator depends on the link variables; its precise 
form will be presented below. 

The average value of the action density, S/V, is directly related to the average plaquette; 
in particular, for the gluonic part we have: 



S — Sg + Sf, 



5G = /3E^G(n), 



(1) 



□ 



with 







2 



As for {Sp/V), it is trivial in any action which is bihncar in the fermion fields [10]; this can 
be seen by rescaling Sp in the fermionic partition function by a factor e : 

Z^{e) = I [D^iV^,] exp (-e^^) = e'^fZ^'ie = 1) (5) 

There follows: 

{T.Ep{x,y)) = -|(^^^)e=i = -^NNf (6) 

We will thus calculate {Eg{0)) in perturbation theory: 

{EG)^cigl + C2g^ + csgl + --- (7) 

The n-loop coefficient can be written as Cn — c^+c^ where is the contribution of diagrams 
without fermion loops and comes from diagrams containing fermions. The coefficients 

have been known for some time up to 3 loops [12,10]; more recently, some estimates 
for higher loops have been produced in 4 dimensions using stochastic methods [13], in the 
absence of fermions. The coefficients are also known to 3 loops for Wilson fermions [10] 
and clover fermions [14]. The task at hand, for overlap fermions, presents many more 
complications, and we will be computing cf^ in the present work. 

The calculation of c„ proceeds most conveniently by computing first the free energy 
— {lnZ)/V, where Z is the full partition function 

Z = j[VUVi^{D^I;,\eM-S). (8) 
The average of Eq is then extracted as follows 

In particular, the perturbative expansion of {lnZ)/V : 

3{N^ -I) , ^ di d 



{\nZ)/V = constant + ( ^-^ ^ ln/3 + ^ + ^ + • • • (10) 



leads immediately to the relations: C2 = di/{2AN'^), C3 = d2/{24N^). 

A total of 8 Feynman diagrams contribute to the present calculation, up to 2-loop order; 
these are shown in Figure 1. 

All the algebra necessary to carry out this calculation was done with our computer 
package written in Mathematica; to this end, certain extensions of the package were needed, 
in order to incorporate the numerical integration of 2-loop expressions involving overlap 
fermions. The first 6 diagrams in Figure 1 give the purely gluonic contributions to the free 
energy: cf , c^, which are well known (see, e.g., [10]). Their contribution is presented here 
for the sake of completeness: 



8A^ ' 

{n' - 1) (o.0051069297 - (11) 
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Fig. 1 Feynman diagrams contributing to the free energy up to 2 loops. Solid, wavy, dashed lines denote 
fermions, gluons and ghosts, respectively. A solid square stands for a contribution from the measure part 

of the action. 

Diagrams 7 and 8 are much more cumbersome to calculate; they lead to the value of , 
whose dependence on N and Nf can be made explicit: 

Co^- -Nfh2, (12) 

where /12 is the quantity we set out to calculate; it depends only on a parameter Mq which 
appears in the overlap action, to which we now turn. 

The Neuberger-Dirac operator Djv appearing in the fermionic part of the action, Eq.(3), 
has the form [15]: 

X^aDw-Mo, (14) 

where is the Wilson-Dirac operator (with the Wilson parameter r set to its standard 
value, r = 1) 
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(v; + V^) - aV;V J , V^^lj{x) = -[U{x,ii)il}{x + aii)-il^{x)], (15) 

a 



Mo is a real parameter, whose value may be chosen at will; in order for this action to describe 
one massless fermion, with no doubling problem, Mq must lie in the range < Mq < 2. 
Nonperturbatively one expects —rric < Mq < 2, where rric < is the critical mass associated 
with the Wilson- Dirac operator. 

Despite the non-ultralocal nature of the square root [X^X]^/^, the operator is nev- 
ertheless local under reasonable assumptions regarding the gauge configuration [16]. In 
particular, the expansion of Djy to all orders in perturbation theory, while very elaborate, 
poses no conceptual obstacle and can be performed starting from an integral representation 
of the square root [17]. 

The expressions for the overlap propagator and vertices are presented in the Appendix. 
It is worth noting that, by the nature of the overlap operator, the fermion vertex with two 
gluons contains two distinct parts: S2 = + S2/3 (see Appendix A); the first of these is 
pointlike, while the second involves integration over an internal four-momentum, as shown in 
Fig. 5. Consequently, diagram 8 splits into parts 8a, 8b; diagram 8b has the same topology 
as diagram 7. 

The numerical integration over loop momenta was carried out for a wide range of vahics 
for < Mq < 2 on finite lattices of size L < 128, with subsequent extrapolation to infinite 
size. The systematic error resulting from the extrapolation is typically small enough to 
leave the first six significant digits intact. A very slight reduction in precision is noticed as 
Mo — > 2, where the pole structure of the fermionic propagator is expected to change. 

Our results, as a function of Mq, are presented in Figure 2 for each fermionic diagram 
separately; their sum is shown in Figure 3 and listed in Table I. It is worth noting that the 
results are very smooth functions of Mq, despite the potential singularities of the integrands 
at Mq = and Mq = 2. 

To facilitate comparison, we provide below the values of individual diagrams for a par- 
ticular value of Mo: 

hl{Mo=l) = -0.02013079(1), /i^"(Mo=l) = -0.0099775548(1), 
hf{Mo=l) = 0.010764939(1) (16) 
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Fig. 2 Numerical results for each diagram involving fcrmions, as a function of the parameter Mq. The 
values denoted by 8a and 8/3 correspond to the two contributions to diagram 8, coming from vertices 

and S2/3, respectively. 
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Fig. 3 The total contribution of fermionic diagrams to /12 , as a function of the parameter Mq • 
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For ease of reference, we list below the total value of {Eg{0)) ior N — 3 and Nf = 0, 2, 3, 
with some typical values of the action parameters {k is the hopping parameter of the Wilson 
action) : 

{Eoia)) = + ■ 0.033910993 {Nf = 0) 

= ^o/3 + 90 ■ 0.0234689 {Nf = 2 Wilson, r = 1, k = 0.156) 

= ^o/3 + 90 ■ 0.0253139 {Nf = 2 overlap, Mq = 1.0) 

= g^/3 + g^ ■ 0.0226752 {Nf = 2 overlap, Mq = 1.5) (17) 

= g^/3 + g^ ■ 0.0182479 {Nf = 3 Wilson, r = 1, k = 0.156) 

= g^/3 + g^ ■ 0.0210154 {Nf = 3 overlap, Mq = 1.0) 

= ^o/3 + 9o ■ 0.0170573 {Nf = 3 overlap, Mq = 1.5) 

The extension of the above calculation in order to encompass Wilson loops of larger size, 
as may be required for computing the quark- ant iquark potential, presents no additional diffi- 
culties. Further two-loop calculations with the overlap action involve fcrmionic vertices with 
3 and 4 gluons; generating these vertices is conceptually straightforward, though technically 
it is quite involved. We are presently working on these extensions, and expect to provide 
results in a future report. 

APPENDIX A 

We present here the expressions for the overlap propagator and for the fermionic vertices 
containing up to two gluons, following Refs. [17,3]. This Appendix serves for completeness, 
and also in order to correct some typographical errors (which had no consequence on previous 
results) . 

Let us first write down the weak coupling expansion of the Wilson-Dirac operator Dw 
This will be useful for constructing the relevant vertices of Dn ■ We write 

M 

X{q,p) ^ a' ^e-^-«+^^-^'^(Dw(a:,Z/) - — S.,y) 

x,y ^ 

= Xo{p){2T,Y5\q - p) + X,{q,p) + ^2(5,^) + 0{gl), (Al) 

i I 1 

where : Xo(p) = - ^ 7^ sin ap^ + - ^(1 - cos ap^) Mq, (A2) 

a ^ a ^ a 

X^{q,p) = go J d^k5{q-p-k)A^{k)V^,^{p + k/2), (A3) 
Vi,,_,{q) = i'^t, cos aq^ + sin aq^,, 

X2{q,p) = I / ^-^^^(^ -p-k,- k,)A,{k,)AMV,Jp + k,/2 + k,l2l (A4) 
y2,^i{q) = -i'y^asmaq^ + acosaq^. 
The Fourier transform of the Neuberger-Dirac operator takes the form 
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-^D^{q,p) = Do(p)(27r)V(g -p) + E(g,p). (A5) 

IVIq 

The factor 1/Mq in the above equation is inconsequential, since it can be absorbed in the 
definition of ip{x). Dq{p) is the tree level inverse propagator: 

2[u{p)^h{p)\ +2' ^^'^ 



where : = - [ ^ sin^ ap^ + [^(1 - cos ap^) - Mq]^ J 



1/2 



1 1 

Kv) = - - cosap^) Mq. 



The function S(g,p) can be expanded in powers of qq as 



(A7) 
(A8) 



aS(g,p) 



u;(g)+u;(p) 
1 

+a;(p) 



+ 



^2(g,p) 

1 



1 



uj{q)uj{p) 
1 



Xo(g)X|(g,p)Xo(p) 



X 



(27r)4 + a;(p) u{q) + + u{p) 

-Xo(g)xl(g,A;)Xi(fc,p) -Xi(g,fc)Xt(A;)Xi(fc,p) -Xi(g,fc)Xj(fc,p)Xo(p) 



a;(g)a'(A;)a;(p) 



+ 0{gt) 



(A9) 



From S(g,p) one can read off the vertices containing up to two gluons. The vertex with one 
gluon, aEi(g,p) corresponds to the first line of Eq.(A9). 



P 



Fig. 4 Pointlike part of the vertex. 
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Fig. 5 Vertex part with intermediate mo- 
mentum integration 



The fermionic vertex with two gluons splits into a pointlike part aL2a (Fig- 4, second line 
of Eq.(A9)) plus a part involving an internal four- momentum k, which must be integrated 
over, aS2/3 (Fig. 5, last 3 lines of Eq.(A9)). 
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TABLE I. Two- loop fermionic con- 
tribution /i2, for various values of Mq, 
< Mo < 2. 



Mo 


h-2 


0.01 


-0.01021118(5 


0.05 


-0.01050736(5 


0.1 


-0.01088833(5 


0.2 


-0.01168252(5 


0.4 


-0.01338658(3 


0.6 


-0.01523365(1 


0.8 


-0.01721907(1 


1. 


-0.01934341(1 


1.2 


-0.02160970(1 


1.4 


-0.02402085(3 


1.6 


-0.02657516(3 


1.8 


-0.02925538(3 


1.9 


-0.03062501(5 


1.95 


-0.03130974(5 


1.99 


-0.0318533(3) 
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